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Abstract 

By extending the lepton sector of standard model to include one sterile neutrino and two sets 
of new Higgs doublets and right-handed neutrinos, denoted by (771, N\, Ns) and (7721^2,-^4), with 
two Z2 symmetries, the puzzles of neutrino masses, matter-antimatter asymmetry and cosmic- 
ray excess observed by Fermi-LAT and PAMELA can be resolved simultaneously. The characters 
of the model are: (a) neutrino masses arise from type-I and radiative seesaw mechanisms; (b) 
leptogenesis leading to baryon asymmetry at the energy scale of 0(1 — lOTeV) could be realized 
through soft Z2 symmetry breaking effects; and (c) the conditions of small couplings for a long-lived 
dark matter could be achieved naturally through loop corrections due to the same soft symmetry 
breaking effects. The candidate for dark matter in leptophilic decays could fit the Fermi-LAT and 
PAMELA data well. 
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One of the strongest motivations for new physics is due to the unsolvable problems in the 
standard model (SM), such as the origin of neutrino masses, matter-antimatter asymmetry 
and dark matter, which are the phenomena observed well, but less understood. Recently, 
three astonishing experiments PAMELA fl, ATIC Q and Fermi-LAT [3| in cosmic-ray 
measurements reheated the issue of the unknown stuff, in which the first one indicates 
anomaly in positron flux ratio while the later two observe excess in electron+positron (e + e _ ) 
flux. Although Fermi-LAT's results don't display the bump as shown at ATIC in the e~e + 
flux, a clear excess at Fermi with high precision can not be denied. Promptly, it is proposed 
that one of possible resources for generating the exotic events in the flux of electrons and 
positrons is ascribed to the dark matter annihilations and/or decays 

PAMELA observes not only on the positron flux ratio but also on the antiproton flux ratio, 
however, the resulting evidences show no significant excess on the latter. Accordingly, if the 
excess observed at Fermi and PAMELA stems from the dark matter, the dark stuff should 
presumably be leptophilic. Furthermore, if the matter-antimatter asymmetry arises from the 
so-called leptogenesis where the baryon asymmetry originates from the lepton asymmetry , 
the issues of neutrino masses, baryon asymmetry and Fermi/PAMELA puzzles intriguingly 
are all related to the lepton sector. Based on this inference, it is interesting to investigate 
how these unsolved problems can be explained within a simple model uniformly. 

Although dark matter annihilations could provide the source for the Fermi/PAMELA 
anomalies, it is inevitable that an enhanced boost factor of a few orders of magnitude, 
such as Sommerfeld enhancement 8|, near-threshold resonances and dark-onium formation 
jsl, has to be introduced. To evade of introducing unnecessary effects, we will focus on 
the mechanism of dark matter decays. Moreover, to get a long-lived dark matter, say 
O(10 26 s), usually we have to fine-tune either the new couplings to be tiny or the new scale 
in intermediated state to be as large as the scale of grand unified theories (GUTs). Hence, 
in order to fulfill our purpose, the subject we face is not only how the model provides a 
new CP violating (CPV) mechanism in the lepton sector so that lepton number asymmetry 
could be converted to the baryon asymmetry by the nonperturbative sphalerons [lol |. but 
also how a long-lived dark matter exists naturally while solving Fermi/PAMELA anomalies. 

To solve the mentioned problems, we extend the SM by including two extra Higgs doublets 
vT = iViiVi) (i=l,2) with hypercharge Y = —1 and five right-handed neutrinos denoted 
by A?!_4 and N. Besides, we introduce two Z 2 discrete symmetries so that rji and A r 1 3 are 
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transformed by 

Vl^-Vl, #1,3 ^"#1,3, (1) 

while r]2 and N 2 ^ follow 

r] 2 -> -r/ 2 , iV 2i4 -> -iV 2 , 4 (2) 

under the first and second Z 2 symmetry, respectively The unmentioned fields in above 
equations denote invariance in each transformation. In addition, rji^ do not have vacuum 
expectation values (VEVs) when the discrete symmetries are exact. Accordingly, the Yukawa 
sector could be found as 

-C Y = L i Y^H£ jR + LjYfHN + L t y taVl N a 
+L iW2 #/3 + ^N T CN 



■i 



+ J2^N^CN k + H.c, (3) 

k=l 

where L T = (v e ,£) L denotes the SU(2) L doublet leptons, H T = (if , if - ) is the SM Higgs 
doublet, Y^, Y^ and yi a (j3) with i,j = 1 — 3 and a = 1,3(2,4) stand for the Yukawa 
couplings. For simplicity, we have chosen Ni as the diagonalized states before the electroweak 
symmetry breaking (EWSB). Since the right-handed neutrino N directly couples to the SM 
Higgs and ordinary leptons, the masses of neutrinos could be induced after the EWSB. 



Therefore, N could lead to non-zero neutrino masses by the type-I seesaw mechanism [11], 
i.e., we have to set 3> tti^. Moreover, because of interactions associated with charged 
leptons being t^N^f, in order to explain the Fermi/PAMELA data, the decaying dark 
matter should be the fermionic particle. Here, we take N 2 as the candidate. Thus, the mass 
relations are required to obey rn Vl{2) > m Nl3(2) and m N , 2 > m Nl . We will see clearly later 
that if m^ A ^> mjy 3 > m m , the role of N 3 is responsible to the leptonic CP asymmetry 
(CPA) in leptogenesis. 

Since is protected by the ^-symmetry and the decay is forbidden by kinematics, so 
far it is still a stable particle. We note that A3 and A4 are allowed kinematically to decay 
through the channels A 3 — > N 2 £ + £~ and A 4 — > £i] 2 , respectively. In order to make the dark 
matter unstable, we ascribe that the origin of the unstable dark matter is from the Z 2 soft 
breaking terms, given by 

-£soft = /WV3 T C#4 + #.C (4) 



Due to the soft breaking effects being associated with heavier right-handed neutrinos, it is 
clear that although Z 2 symmetries have been broken, the breaking effects will not open a 
sizable decaying channel for N 2 . Hence, the unstable dark matter could be a long-lived one. 

We now start to discuss in turn how the puzzles of small neutrino masses, baryon asym- 
metry and Fermi-LAT/PAMELA are solved in this model. 

Neutrino Masses: By the first term of Eq. ([3]), the neutrino mass matrix could be ex- 
pressed by 

( m ^)J yPC_I ~ -m Di m- N l m ] Djl m Dj = ~j=Y* , (5) 

where we have made the expansion in terms of mo/mN and used (H) = v/V2 as the 
vacuum expectation value (VEV) of the Higgs field. To explain the neutrino masses, with 
ni) ~ 0(1) we see that ~ 3 x 10 15 GeV/(m l/ /leV) as expected by the conventional 
seesaw model. Beside the type-I seesaw mechanism, the masses of neutrinos could be also 
induced by radiative corrections in which the corresponding Feynman diagram is illustrated 
in Fig.[3](a). With the relevant quartic term in the Higgs potential, written by A5/2 (H^rji^) 2 , 
the mass matrix through one-loop is given by 1^] 

/ N rad A 5 w 2 / yayjim Nl 



,,J 8tt 2 V m 2 



yi3yj3m m yi2yj2m N2 \ 



(6) 



where we have adopted m Vl , 2) > Ttjv 13(2) - Because iV 4 is much heavier than other particles, 
we have ignored its contributions. With m Vl ~ 8 TeV, m m ~ ~ 4 TeV, ~ mjv 2 ~ 
O(lTeV), A5 ~ O(f0 -4 ) and y^ ~ O(f0~ 2 ), the neutrino mass by the dominant effects could 
be O(leV). Intriguingly, if we combine the type-I and radiative seesaw mechanisms together, 
it is found that some cancelations could occur between both. As a result, the constraint on 
the couplings of A 5 and ya(2,3) could be somewhat relaxed. Besides, the m u ~ leV resulted 
in each seesaw mechanism can be reduced to 0(0.01 — O.leV) by the cancelations. 

Leptogenesis: Since iV 3 plays the role of the lepton asymmetry at the energy scale of 
0(1 — lOTeV), to satisfy the out-of-equilibrium condition, the decay rate of A3 should be 
less than the Hubble constant H at the temperature of mjv 3 , given by 13[] 



rjv 3 < H(T = m N3 ) = \/ -— | r=ro (7) 

40 Mpi anck 
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where g* ~ 100 denotes the number of active degrees of freedom and Mpi anc k ~ 10 19 GeV 
is the Planck scale. As mentioned early, N3 could decay via the three-body channel of 
— > Ni£ + £~. We check that by the phase space suppression and with y^yn ~ 10~ 5 , the 
rate for the three-body decay is around a factor of 5 smaller than Hubble constant. In 
addition, by the soft breaking effects of Eq. (j3J), N$ will decay to irji through the mixing 
with N4, sketched by Fig. [U(b). The rate will depend on the parameters of /134, rriN 4 and y^. 
By calculating the two-body decay rate of Fig. [0(b) and with rriN 3 = 4 TeV, the constraint 
on the parameters is 



/^34 

-yi4 



< 5 x lO^ 7 . (8) 

By employing small ^34/772^4, the decay rate less than the Hubble constant for A3 of O(TeV) 
can be naturally accomplished. 

Another necessary condition to achieve baryon asymmetry is the CPA, defined by 

A = T(N 3 -> £4) - T(N 3 -> i%d 

A nonzero direct CPA should involve CPV and CP conserving (CPC) phases simultaneously. 
Here, the complex Yukawa couplings yi a can provide the new CPV source. Therefore, we 
need a new mechanism to generate the physical CPC phase. It has been well known that 
one-loop effects could produce the CPC phase when the on-shell condition of particle in the 
internal loop is satisfied. In the model the CPC phases could be induced via self-energy 
[bll ] and vertex corrections {3] illustrated in Fig. [TJc) and (d). Therefore, the CPAs via 




FIG. 1: Feynman diagrams of (a) Neutrino masses by radiative corrections and (b), (c) and (d) 
the N 3 — > ir\2 decay at tree and loop levels. 
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self-energy (S) and vertex (V) diagrams are respectively found to be [1 



A CP (N 3 ) « sin 5 



1 • 



x 



x 2 [(1 + x 2 ) log(l + l/x 2 ) - 1] 



(10) 



with %2 = (m7v 2 /m7v 3 ) 2 and r y sin<5 = Irn{y* i y i 2) 2 /\yi^\ 2 . 5 is taken as the CPV phase and 
the repeating index % stands for the summation in lepton flavors. Followin g th e relation of 



the observed baryon asymmetry and the lepton asymmetry, formulated by 

fn B \ 2Sn B — n L n L n. 



16| 



-A 



CP 



(11) 



V s Jobs 79 s s 2s 

where s and ?t<b(l, 7 ) are the entropy and baryon (lepton,7) density of the universe, respec 
tively, we obtain 



S ' obs 



I2g, 



(12) 



where we have used s = g*T 3 (27i 2 /45) and 



2T 3 /vr 2 . With | sin <5 1 = 0.3, the baryon 



asymmetry as a function of ry is presented in Fig. [2], where the solid, dashed, dotted and 
dash-dotted lines represent N 2 = 2, 2.5, 3.0 and 3.5 TeV, respectively, with N 3 = 4 TeV. 
According to the results, ry ~ 1 2/*2 1 ~ O(10~ 3 ) not only satisfies the criterion for radiative 




FIG. 2: Baryon asymmetry by leptogenesis, where the solid, dashed, dotted and dash-dotted lines 
represent N 2 = 2, 2.5, 3.0 and 3.5 TeV, respectively, with N$ = 4 TeV. 

neutrino masses but also explains the baryon asymmetry in leptogenesis. 

Fermi-L AT/PAMELA: Due to being heavier than other particles, the soft break- 
ing interactions introduced in Eq. (jlj) cannot directly lead to N 2 decays. However, by 
combining the soft breaking effects with Yukawa couplings, we find that a finite dimension-4 
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hard operator for the mixture of rji and 7/2 could be induced at one-loop level, sketched in 
Fig. [31 The resultant contributions to the scalar potential is given by 

H r,\ 
\ Li / 

X 

/ ^ \ . 

FIG. 3: Box diagram for dimension-4 term (i]\h) (ij/^Hj. 



N 
N 



SV haxd = C h (vlH^j (jilH) + H.c, (13) 

with 

(14) 

where we have neglected the terms associated with mN 3i /mN and used the Einstein sum- 
mation convention for the indices % and j. If we simply ignore the summation in indices % 
and j, by combining with Eq. (jHJ), the mixing effect of r]i and 772, denoted by //f 2 > could be 
simplified to be 

Ml2 ~ T^( m u)ijy*3y^^34 ■ (15) 
107T 

With using y i3 ~ 10~ 2 and mjv 4 = 5 x 10 4 TeV and following Eqs. (jSJ) and (|T5|) . we immedi- 
ately have /ia 2 2 < 1-5 x 10 _12 GeV 2 . Since Z2 symmetries will be restored while /134 vanishes, 
the smallness of /x 2 2 is still nature in technique. Consequently, the lifetime for dictated 
by iV2 — > N\ll can be found by 

TN ~ (16) 

d|S/»22/,-i| m^ 2 (1 - m 2 N Jm 2 N2 ) 

Adopting m Vl , 2) = 8(4) TeV, rriN 1(2) = 0.5(2) TeV, yaSta ~ 7.5 x 10~ 6 and the obtained 
upper value of /i 2 2 , we get r^ 2 « 2.5 x 10 26 s. Clearly, the values of the parameters in 
the model could be compatible with the phenomena of neutrino masses, leptogenesis and 
long-lived dark matter. 



Since the long-lived dark matter could only decay in leptophilic, for illustration, we simply 
consider the electron-positron pair and one electron (positron) and antimuon (muon) as the 
possible final states, in which antimuon (muon) will further decay to positron (electron) 
by SM weak interactions. Due to three possible channels involved, where the associated 
parameters are proportional to |2/ e 22/ei 1 2 ? 1 2/e22/^i 1 2 and |y M 2?/eil 2 ; ^he energy spectrum for 
electron (positron) can be formulated by 

1 



dN e ± 
dE 



2 dY ee 2 dY e ^ 2 fie 

N \ Ca ltE +C ?llE + °^ltE 



(17) 



with 



dE 

d^ e/i 

~Je 

d^ fie 

~dE 



[ ^Am 2 21 - ^m N2 E 



Fe,{E) 
Ffie(E) 



m N „/2 



dE^F^E^) 
dE^Fe^Efj) 



d\Mtf{E„E) 
dE 



d\Mtf{E„E) 
dE 



where N is the normalization, c 2 a + c| + c 2 = 1 , Am?, 



21 



m N 2 m N 1 1 



F e ,(E) = E 2 (A 



rn 



21 



2m N2 E), F^E) = E 2 (Am 2 21 /2 - 2/3m N2 E) and 



dlM^jE^E) = GlJ_U 
67r 3 I 3 



dE 



E d (E,-\p,\r-[^ 



in 



2 \ 3 



3 2 
--7717. 

2 



e 2 (e, - \p,\r - [ --f 



m, 



(18) 



For dT efl /dE and dT^/dE, the allowed energy range for muon and electron are found to 
be m M < En < m x /2 and < E < m 2 l /2(E fJi — respectively. Since the experiments 

measure the flux of cosmic-rays, the used formalism to estimate the flux from the new source 
is given by 

m A r 2 /2 

(19) 



1 



47T mN 2 T N 2 



dE'G(E,E') dNe± 



dE' 



with_c being the speed of light. For numerical estimations, we adopt the result parametrized 
by 



in 



G{E,E')~^-exp[a + b{E 
E z 



5-1 



E' 5 ~ 1 )]8(E' - E) [cm~ 3 s] 



(20) 
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with a 



-0.9809, b = -1.1456 and 5 = 0.46. 



For including the primary and secondary electrons and secondary positrons, we use the 
parametrizations, given by 



<f> s ^(E) 



0.16E 



-i.i 



1 + HE - 9 + 3.2E 2 - 15 
0.7E - 7 



1 + 110E 1 - 5 + 600£ 2 - 9 + 580E 4 - 2 
A.5E - 7 



[GeV-^m^s- 1 ^- 1 ], 



[GeVW^s^sr- 1 ] 



(21) 



1 + 650£ 2 - 3 + 1500£ 4 - 2 

where $P rim ( sec ) denotes the primary (secondary) cosmic ray. Accordingly, the total electron 
and positron fluxes are defined by 



$ pnm $S ec 
e e 



/V 2 



N 2 



(22) 



Here, according to Refs. [19| and [20J, we have regarded the normalization of the primary 
electron flux to be undetermined and parametrized by the parameter of k. The value of k is 
chosen to fit the data. Before introducing the source of the primary positron, k is set to be 
0.8. Taking rriN 2(1) = 2(0.2) TeV and k = 0.65, we present the e~e + flux and ratio of fluxes 
e + / (e~ + e + ) by N 2 decays in Fig. HI in which the thick solid, dashed and dash-dotted lines 



stand for (c 2 , c|, c 2 ) 



'1,0,0), (0,1,0) and (0,0,1), respectively. Obviously, with proper 



" 1000 



s 

u 

o 



"fa 




10 



* Fermi-LAT 



- Background 





J I L 



PAMELA 

Background 

J 



10 100 1000 

E[GeV] 



10 



100 

E[GeV] 



1000 



FIG. 4: electron +positron flux (left) and the ratio of positron flux to electron +positron flux (right), 
where the thick solid, dashed and dash-dotted lines represent (c 2 ,Co,c 2 ) = (1,0,0), (0,1,0) and 



(0, 0, 1), respectively. 



values for the parameters, the decaying dark matter could fit the Fermi-LAT and PAMELA 
data well simultaneously. 
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Inspired by the recent data measured by Fermi-LAT and PAMELA, we have found that 
the puzzles of small neutrino masses and baryon asymmetry as well as dark matter can 
be solved uniformly by extending the lepton sector of the SM with two Z 2 symmetries. In 
the proposed model, the conventional type-I and radiative corrections seesaw mechanisms 
coexist to generate the masses of neutrinos. The matter-antimatter asymmetry originated 
from leptogenesis could be accomplished by the introduced new stuff associated with soft 
Z2 breaking terms. With the same soft breaking effects, the condition of small couplings for 
a long-lived decaying dark matter can be realized naturally through radiative corrections. 
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